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We have studied quantum wires using the Green's function technique and the density-functional 
theory, calculating the electronic structure and the conductance. All the numerics are implemented 
using the finite-element method with a high-order polynomial basis. For short wires, i.e. quantum 
point contacts, the zero-bias conductance shows, as a function of the gate voltage and at a finite 
temperature, a plateau at around 0.7Go. (Go = 2e^ /h is the quantum conductance). The behavior, 
which is caused in our mean-field model by spontaneous spin polarization in the constriction, is 
reminiscent of the so-called 0.7-anomaly observed in experiments. In our model the temperature 
and the wire length affect the conductance-gate voltage curves in the same way as in the measured 
data. 
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I. INTRODUCTION 

Two-dimensional (2D) nanostructures can be 
fabricated at semiconductor interfaces (such as 
GaAs/AlGaAs), using lithographic techniques and gate 
electrodes [1-3] . In these structures conducting electrons 
form a quasi-2D electron gas at the interface. Doping 
impurities are separated from the 2D electron gas mak- 
ing the transport baUistic over the nanostructure. This 
enhances the importance of quantum mechanical effects. 

An interesting 2D nanostructure is the quantum wire 
(QW), which is a laterally narrow electron pathway con- 
necting two (infinite) electrodes. Because of the ballistic 
electron transport the conductance is quantized in units 
of the quantum Gq — 2e^/h. The 2D electron states are 
quantized in the direction perpendicular to the axis of 
the wire. Each occupied perpendicular state gives rise 
to one conducting mode adding one quantum to the con- 
ductance. In a typical measurement, the conductance is 
recorded as a function of the gate voltage which lowers 
the potential in the wire region, influencing its width or 
the average electron density. The increase in the number 
of conducting modes results in a staircase as a function 
of the gate voltage. Recently, Kane et al. [1] and Reilly 
et al. [4] have measured conductances of QW's with dif- 
ferent geometries. The quantization of the conductance 
is clearly seen in the measurements. 

A short quantum wire can be characterized as a quan- 
tum point contact (QPC). QPC's exhibit the so-called 
0.7-anomaly which is a small plateau in conductance 
around 0.7-0.5 Gq as a function of the gate voltage [2]. 
The plateau moves from 0.7 Go to 0.5 Gq as the electron 
density in the quantum wire is lowered or the QPC is 
in a strong magnetic field. The plateau becomes more 
pronounced when temperature increases from the milli- 
Kelvin region to a few Kelvin. At higher temperatures 
the anomaly disappears. The length of the wire affects 
also the anomaly [4] so that in long wires the anomaly is 



around 0.5 Go and in shorter ones around 0.7 Go. The 
measurements show a short plateau also within the next 
quantum step around 1.5 ... 1.76 Go. 

There are several different explanations for the 0.7 
anomaly. One of them [5,6] is based on the Kondo model 
where an unpaired electron is localized at the QPC. In the 
ground state this electron is in the spin-singlet state with 
the scattering electrons from the leads. The spin coupling 
results in a high density of states, a Kondo resonance, at 
the Fermi level and thus in an enhanced conductance 
around zero bias voltage. This zero-bias anomaly is sim- 
ilar to the Kondo effect in single-electron transistors [7] . 
In the case of the QPC the Kondo model seemingly ex- 
plains the plateaus and their behavior as a function of 
the temperature and the external magnetic field [5,6] 

Another explanation is provided by the semi-empirical 
model by Reilly et al. [4]. In this model the plateau is 
caused by the local spin polarization of the electron gas 
at the QPC. For small electron densities at the QPC the 
polarization vanishes, but when the electron density is 
increased using the gate voltage a spin gap opens and 
causes the polarization. The opening of the spin-gap de- 
pends on the QPC configuration and temperature, and 
in different cases the conductance plateaus appear at dif- 
ferent values. A spin-polarization model has also been 
presented by Berggren and Yakimenko [8] . They used the 
density-functional theory (DFT) to calculate the electron 
density in a quantum wire as function of the gate volt- 
age. The conductance was calculated using the Biittiker 
model [9] with parameters obtained by fitting the effec- 
tive DFT potential in the middle of the QW. Also Meir 
et al. [6] and Hirosc et al. [5] reported DFT calculations 
for a QW in order to model the QPC resonance states 
needed in the Kondo model. Also recent measurements 
utilising the resonant interaction between couplet QW's 
show evidence of the formation of localized magnetic mo- 
ment in a constriction of the 2D electron gas [10,11]. 

In this work we investigate the spontaneous spin- 



polarization model for QW's. Our aim is to make within 
this model accurate and realistic calculations for the con- 
ductance, in order to study how far the model can de- 
scribe the experimental findings. We use the DFT and 
the Green's functions techniques (see, for example Xue 
et al. [12]) to calculate the electron density and the con- 
ductance. For the electron spin densities the use of the 
Green's function technique is computationally more de- 
manding than solving for the wave functions of a finite 
system. However, it has important advantages. First, an 
infinite system without any artificial periodicity can be 
treated. It is also possible to add the bias voltage between 
the electrodes and calculate the charge density and the 
current through the structure self-consistently. We solve 
for the DFT-Green's function equations numerically by 
the finite-element method using our recent implementa- 
tion [13]. 

We use the effective atomic units which are derived by 
setting the fundamental constants e = h = rrie = 1, and 
the material constants m* = € — 1. m* and e are the 
relative effective electron mass and the relative dielec- 
tric constant respectively. The effective atomic units are 
transformed to the usual atomic and SI units with the 
relations 



Length: la^ = l_l-ao 
Energy: 1 Ha* = l^f Ha 
Current: la.u.* — l^a.u. 



10.0307 nm 
11.3079 mcV 
2.7512 ^A. 



Above, the numbers on the right-hand side are obtained 
by using the parameters for GaAs, i.e. m* = 0.067 and 
e = 12.7. They are used in illustrating our results. Ha 
denotes the Hartree energy unit. 



of a gate voltage. The denser the background charge the 
deeper the effective potential is, i.e. 
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(1) 



where r — {rx,ry} are the 2D coordinates. We calculate 
Vg at the midpoint of the wire and use this value as the 
gate voltage in our illustrations below. The dependence 
of the gate voltage can be written in the form 



Vn^a 



JPJ^ 



(2) 



where the coefficient Cj depends on the shape of the wire. 
Cj = 218mVa;5^ 23AmVaf and 248mVa*2 for wires with 
the width W of Soq and lengths L of Ga^, la^ and Sap, 
respectively. In chapter III we denote these wires as D, 
E and F. 




FIG. 1. 2D quantum wire between two electrodes. The 
gray areas denote the rigid positive background charge. The 
electrodes continue to the infinity outside the calculation re- 
gion. The uniform background charge density is varied in the 
dark-gray wire region in order to model the effects of the gate 
voltage. 



II. MODEL AND THEORETICAL METHODS 

Our strictly 2D model for a QW is shown in Fig. 1. 
We use a simplified geometry in order to see the effects 
due to different shapes of the wire. The calculation area, 
O, consists of the QW and some parts of the electrodes. 

The semi-infinite electrodes consist of the posi- 
tive background charge with the constant density of 
0.2(aQ)^^ « 2 X lO^^e/cm^ and the neutralizing 2D elec- 
tron gas with the density p{r). 

At both sides of the electrodes and the wire we include 
some empty vacuum, and the electron density is required 
to vanish at the outer boundaries parallel to the wires. 

The electron density at the axis of the electrode is 
then typical for GaAs/GaAlAs systems and the Fermi 
energy (the width of the occupied energy band) is 
0.63 Ha=7.1 meV. The neutrality makes the calculation 
of the electrostatic potential in the self-consistent DFT 
calculations possible. The QW is also modeled using a 
rigid uniform positive background charge. We control 
the positive charge density pj in the wire region ^qw 
(the darker area in Fig. 1) in order to mimic the effects 



We use the DFT within the local density- 
approximation to calculate the electron density. In the 
self-consistent iterations the electron density is first cal- 
culated using Green's functions. The electron density is 
then used to calculate the effective potential Fg// as 



Veffir) 



p{r')-p+{r') 
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dr' + V^cir) + Vm. 
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Above, the first term on the right-hand side is the 
Coulomb potential, computed as an integral and not from 
the Poisson equation because of the 2D calculation space. 
The electron density in the electrodes ^l/r. is assumed 
to coincide with that in an infinite uniform wire. Then 
we can include the outside regions into the Coulomb inte- 
gral within the buffer regions fi^ and fl'j^ which are large 
enough. For the exchange-correlation potential, Vxc, we 
use the 2D-LDA functional by Attaccalite et al. [14,15]. 
The potential Vbias is a linear ramp taking care of the 
boundary conditions under a finite bias so that the po- 
tential in the left electrode is shifted relative to that in 
the right electrode by a given bias voltage Vsd ■ The effec- 
tive potential determines through the Green's functions a 



new electron density and the procedure is repeated until 
convergence. 

The retarded Green's function is obtained by solving 
the equation 

(c + iv2 - V,ff{r))G^{ry-uj) ^ S{r r'), (4) 

where to is the electron energy. At the boundaries dQ^/n 
we use open boundary conditions given by the Dirichlet- 
to-Neumann map on the boundary. This means that 
there are no reflections. One can next calculate the so- 
called lesser Green's function, G'^(r, r'; w). For zero bias, 
when the Fermi functions in the left and right electrodes 
are identical, i.e. fii^) = Ir{'-^): G^ is obtained as 



G<(r,r';^) = 2/i/«HG'-(r,r';^). 



(5) 



In the non-equilibrium situation when fhi^) 7^ fni'-^) wc 
have 

G<{r,r';uj) = 
- ifRi^) 



G''{r,rR]uj)VR{rR,r'j^;uj) 
xG''{r'R,r';Lu)drndr'j, 






(6) 



Above, Ti^/n are functions connecting the calculation 
area to the outside electrodes. They have nonzero values 
only on the boundaries dQ.i^m. The imaginary part of 
G< is related to the density of states, and the electron 
density is calculated by integrating over the energy 



-1 /■°° 
p{r) = / Im(G^(r, r: a;))da 



(7) 



The local density of states (LDOS) in the QW region is 
calculated as 



9Qw{^)^^r-j Im(G<(r,r;w))d 



(8) 



LDOS is a continuous function in our calculations and 
thus the finite-size effects are small. 

The electron tunneling probability T{lo) is calculated 
as a function of the electron energy using the Green's 
functions as 

TH=I III ri(ri,r^;c^)G'-K,r^;^) 
JanL Joul Jdnn. Jdnn 

xTR{rR,r'i^;u)G''{r'n,rL;io)drLdr'LdrRdr'ji, (9) 

Where G" is Green's advanced function. Thus one 
needs the Green's function values only at the boundaries 



as 



L/R- 



For a finite bias the electric current is calculated 



/ = - / T{cu) ihiuj) - fR{uj)) duj, (10) 

n 



where the Fermi functions /^ and fR are shifted with re- 
spect to each other by the bias voltage Vsd ■ In the zero 
bias limit fL{^) = ./i?(^) — /(^) ^nd one obtains the 
linear-response conductance 



G 



1 iuj) ; diV. 

dcj 



(11) 



At zero temperature the conductance is simply T(a;/), 
where cof is the Fermi energy. At a finite temperature 
also electron states with energies near the Fermi level 
contribute to the conductance, as the derivate of the 
Fermi function /(w) differs from the delta function. A fi- 
nite temperature infiuences the solution also through the 
electron density (cf. Eqs (5) and (6)). Below we show 
also differential conductances corresponding to given bias 
voltages. At zero temperature they can be computed 
from the approximation 



dl 

dv7 



sd 



I{Vsd + 5V)-I{Vsd) 
6V 



(12) 



where 5V is a small increment in the bias voltage. 

We have numerically implemented the non-equilibrium 
DFT-based scheme using the finite-element method as 
explained earlier [13]. The size of the matrices to be in- 
verted is determined by the number of the finite-element 
basis functions needed in the calculations. We have im- 
plemented 2D high-order polynomial bases [16] up to the 
fourth order in order to reduce the basis size. In a typi- 
cal calculation the number of basis functions ranges from 
2800 for to high-order polynomials to 5500 for low-order 
polynomials. 



III. RESULTS AND DISCUSSION 

A. Spin polarization in the QW 

In this model, there are electronic resonance states in 
the QW, which are reflected as peaks in the LDOS (see 
Fig. 7 below). The peaks are broader in short and wide 
QW's than in long and thin ones, because short QW's 
are more strongly connected to the electrodes. If the res- 
onance peaks are narrow enough, a spontaneous spin po- 
larization may occur with a spin gap opening between the 
spin-up and spin-down states below and above the Fermi 
level, respectively. The polarization occurs in a limited 
range of gate voltages. At small gate voltages there is 
no polarization and it disappears again at large voltages. 
This is in accordance with the density-dependent opening 
of the spin gap speculated by Reilly et al. [4]. A solu- 
tion with spin polarization in the QW is shown in Fig. 2 
which gives the total electron density and the difference 
between the spin-up and spin-down densities. The elec- 
tron wave functions have in this case no nodes in the 
direction perpendicular to the QW, so that there is only 
a single conducting mode. 




FIG. 2. (a) Total electron density and (b) difference be- 
tween the spin-up and spin down electron densities at the 
zero temperature for the system with the dimensions S=5 aj, 
L=7 a^, W=20 aj, A=47 aj (see Fig. 1) and for the gate 
voltage Vg^ 12 A mV 

Berggren and Yakimenko [8] have studied the spin po- 
larization in QW's. In their calculations the electrodes 
are modeled by two large quantum dots which are con- 
nected by the QW. The potential and the electron density 
inside the QW are controlled by a gate voltage. The elec- 
tron density is calculated from the wave functions using 
the DFT within the LDA. They obtain a spin polariza- 
tion which closely resembles our results. The spread of 
the polarization in the electrodes is somewhat larger than 
in our cases which show only a small oscillating polariza- 
tion outside the QW. 

Berggren and Yakimenko do not report on the details 
of the resonance peaks in the DOS. The resonance state 
parameters are crucial in the Kondo model. Therefore 
Meir et al. [6] have used the DFT to model the QPC. In 
their calculations the semi-infinite electrodes are mod- 
eled by wires with a parabolic confinement potential per- 
pendicular to the wire axis, and a raise in the external 
potential forms the QPC as a saddle-point. Due to the 
different potential construction the resonance peaks ob- 
tained by Meir et al. are remarkably wider than those 
in our calculations. A similar QPC construction is used 
also in the calculations by Hirose et al. [18]. The polar- 
izations they obtain are similar to our results for long 
QW's. 



B. Conductance as function of the width of the wire 

The conductances of three QW's A, B and C with dif- 
ferent widths are shown in Fig. 3 as a function of the gate 
voltage at zero temperature. The QW's have the same 
length L=7 ag but for A, B, and C the widths S=5 a|5, 
6 a^ and 10 ag, respectively. The staircase quantizations 
of the conductances are clearly seen. For the widest wire 
C the perpendicular states are denser in energy than for 
the thinner wires A and B, so that the conductance steps 



are shortest for C. The first plateau for C is also dis- 
turbed by the states decaying from the electrodes to the 
QW. Wires A and B are connected more weakly to the 
electrodes and for them the first plateau is clearer. A 
weak connection causes also the spin-polarization in wire 
A, seen as a short plateau clearly below 1 Gg. Wire C 
is not polarized and in wire B the polarization is not 
strong enough to cause plateaus deviating from integer 
of Gg values. 
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FIG. 3. Conductance at zero temperature as a function of 
the gate voltage for three wires with the length L=7 a2 and 
widths S=5 a5 (A), 6 ao (B) and 10 aS (C). The width of the 
electrodes is W=20 aS, and the length of the computational 
area is A=47 ag (see Fig. 1.). 

The gate voltage drives the electrostatic potential quite 
uniformly within the QW. In experiments, side gates are 
often used to control the potential [1]. These gates not 
only increase the potential in the QW, but they also make 
it narrower. The result is a conductance staircase as a 
function of the gate voltage. The first measured plateaus 
corresponding to IGg are pronounced [1] in contrast to 
our shoulder- like first plateau for wire C. Below we are 
concerned mainly with the conductance plateaus below 
IGg and therefore we discuss wires with the width S=5 a^ 
(wire A in Fig. 3). 



C. Conductance as function of the length of the wire 

The conductances of the wires with the width S=5 a^ 
and with different lengths L arc shown in Fig. 4 as a 
function of the gate voltage and at zero temperature. 
The figure shows clearly the effect of the electrode-wire 
connection. The long wires have clear peaks due to the 
resonances. The heights of the peaks are 0.5 Gg, mean- 
ing that only a single electron polarized mode contributes 
to them. The wires with the lengths L=6...8 ag exhibit 
resonances which are just narrow enough for the spin 
polarization to appear. The length dependence of the 



conductance among these three wires is in quahtative 
agreement with the recent measurements for QW's by 
ReiUy et al. [4], although the wires in the experiments 
are clearly longer than in our calculations. 



the polarization reaches its maximal value. 




D. Temperature dependence of the conductance 



10 20 

Gate voltage (mV) 

FIG. 4. Conductance as a function of the gate voltage for 
QW's with the width S=5 ao and with different lengths L. 
The width of the electrodes is W=20 ao, and the length of 
the computational area is A=47 aj (see also Fig. 1). The 
successive curves have been shifted by 0.5 Go- The conduc- 
tance of the wire with L=10 ao is decomposed into spin-up 
and spin-down contributions (dotted lines). 

Berggrcn and Yakimcnko [8] have also calculated the 
conductance of their QW system using a rough scheme. 
They fit parabolic curves to the effective potential in the 
middle of the QW. The fitting parameters are then used 
in the Biittiker model for the conductance. The results 
have the same overall features as ours for the longest 
wires in Fig. 4, but the changes in the polarization and 
in the conductance are less abrupt in our calculations. 
Berggren and Yakimenko argue that the rapid changes 
are due to the finite-size effects. Therefore our results 
describing infinite systems seem to support their state- 
ment. The conductance-gate voltage curves by Berggren 
and Yakimcnko show two plateaus, one plateau when the 
polarization appears and another one at ~0.75 Go, just 
before the polarization disappears. This result resembles 
our curves for the long wire with L=10 Oq. However, 
there is a notable difference that in our calculations for 
the long wires the polarization at low gate voltages is first 
nearly perfect whereas that by Berggrcn and Yakimcnko 
is first absent and appears then as a partial polarization. 
In our conductance curves only one plateau is seen when 



The effect of the temperature on the conductance be- 
havior of three QPC-like wires D, E and F is shown in 
Fig 5. At zero temperature, wire D (L=6 aj^) shows 
no plateau, whereas wire E (L=7 ag) has a plateau at 



-0.7 G5 and whe F (L= 



=8aS, 



at ~0.5 Go. When tem- 



perature increases the plateaus below IGo in wires E 
and F shift down wards and become smoother. Wire D 
shows a weak temperature dependence so that the slope 
at '-~-0.7Go decreases. 

Reilly et al. [19] have presented a phcnomcnological 
model for the appearance and the location of the con- 
ductance plateaus. They report on the measurements 
of three QPC's and explain their conductance behaviors. 
Two of the QPC's have plateaus at ~0.5 Go and ^0.7 Go 
and the third one has a very weak plateau at ^^0.6 Go. 
The QPC's, D, E and F in Fig. 5 show similar gate- 
voltage and temperature dependences. 

In the model of Reilly et al. the relevant parameter 
is the ratio between the spin gap width and the thermal 
energy ksT. According to our calculations the temper- 
ature broadening in the Fermi functions is small when 
compared to the resonance peak widths arising from the 
interactions with the electrodes. Because of this, ksT 
has to be replaced by the resonance peak width as the 
relevant parameter. Then the model can be used to ex- 
plain the conductance- gate voltage curves also at zero 
temperature. 
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FIG. 5. Conductance as a function of the gate voltage for 
QW's with the width S=5 aj and lengths L=6 aj (D), 7 a^ 
(E) and 8 aj) (F) at temperatures of K (solid curve), 2 K 
(dashed curve) and 4 K (dotted curve). 



FIG. 6. Conductances due to spin-up (dashed curves) and 
spin-down (solid-curves) electrons as a function of the gate 
voltage and at temperatures of K and 2 K. The results are 
for wire E (Fig. 5) with the length L=6.5 Og. 



Reilly et al. [19] have given three scenarios based on 
their model. In their scenario I the spin polarization in- 
creases rapidly with the increase of the electron density 
at the QPC. Then a plateau appears at ~0.5 Gq. This 
corresponds to our wire F, in which the resonance peaks 
are narrowest. In the scenario II by Reilly et al. the tem- 
perature broadcnings (the widths of the resonance peaks 
in our generalization) are comparable to the spin gap. 
This corresponds to our wire E. When the spin gap be- 
gins to open at a certain gate voltage also the spin-down 
resonance states are partially populated. This situation 
remains when the density increases and results in the con- 
ductance plateau at ^0.7 Gq. In scenario III by Reilly et 
al. the spin splitting is weak in comparison to the width 
of the resonance peaks. This corresponds to our wire D 
for which no plateau is seen at low temperatures. 



The reason for the different temperature behaviors of 
the QW's in our model can be seen from the curves in 
which the conductances caused by spin-up and spin-down 
electrons are separated. For wire E this is shown in Fig. 6. 
As temperature increases from K to 2 K spin polariza- 
tion increases for gate voltages around the middle of the 
plateau below IGq. At the same time, the resonance 
peaks also become wider because more states contribute 
to the conductance (see Eq. (11)). The same behavior is 
seen also for wires D and F. The reason for the increase 
in the polarization is seen in the LDOS for wire E in 
Fig. 7. When the temperature rises the electron density 
increases in the QW due to the resonances near the Fermi 
level. Then the decrease in the exchange-correlation en- 
ergy opens the spin gap, as can be seen in Fig. 7, and 
the polarization increases. The effect is clearer for wire 
E than for wires D and F. In wire D the electron density 
does not increase as fast as in wire E because the reso- 
nance peaks are wider. Wire F has a strong polarization 
already at zero temperature and therefore it cannot show 
an increase as large as wire E. 
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FIG. 7. LDOS corresponding to the QW region (Eq. 8) 
wire E with the width S=5 ag and length L=7 Og. The de- 
composition to the spin-up and spin-down states at the gate 
voltage of 14 mV and at the temperatures of K (solid-curve) 
and 2 K (dashed-curve) are given. The dotted line denotes 
the Fermi-level. 



The differential conductance calculated using Eq. (11) 
is shown in Fig. 8 as a function of the gate voltage. The 
increase in the applied bias voltage Vsd mainly increases 
the conductance but at the same time it also curtails the 
conductance plateaus. This is not exactly in agreement 
with the measurements which show that the conductance 
plateau below 1 Go rises with the increasing bias volt- 
age [20,5]. According to our calculations the raise of Vsd 
causes also the decrease in the spin polarization in the 
QW. 

The measured differential conductance suppresses 
rapidly as a function of the source-drain voltage. This 
is called as the zero-bias anomaly. Our mean-field model 
cannot produce this kind of behavior whereas the Kondo 
model [6] gives an explanation for the zero-bias anomaly. 
In the Kondo model an unpaired electron couples to the 
electrons of both leads, resulting in an increased LDOS 
(Kondo resonance) at the Fermi levels of the electrodes 
and in an enhanced conductance. When the separation 
between the Fermi levels increases with bias voltage, elec- 
trons can no longer resonantly tunnel through the QPC 
which leading to the suppression in the conductance. 



Meir tt al. [6] have modeled a QW's using the Kondo 
model [6] as explained in the Introduction. In their model 
the Anderson Hamiltonian with parameters controlling 
the properties of the QPC is used. In their example the 
conductance plateau is smooth at low temperatures and 
located around ^0.9 Gq. As the temperature increases 
the plateau becomes wider and moves down to '-^^0.7 Gq. 
The behavior is similar to that observed here for wire E. 



E. Effect of the bias voltage 
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FIG. 8. Differential conductance of wire E with the width 
S=5 ao and length L=7ao as a function of the gate voltage 
and at zero temperature. Results for the bias voltages of mV 
(dotted curve), 0.23 mV (dashed curve) and 0.46 mV (solid 
curve) are given. 



IV. CONCLUSIONS 



We have used the density-functional theory and the 
Green's function formalism to model electronic trans- 
port through ballistic 2D quantum wires. Our results 
show the appearance of a spontaneous spin polarization 
in the wire. Under proper conditions, the spin polariza- 
tion causes a plateau at around 0.7 Gq in the conduc- 
tance as a function of the gate voltage. The calculations 
explain measurements by Kane et al. [1] for the move- 
ment of the conductance plateau as function of the wire 
length. Also the temperature dependence is qualitatively 
similar to the experimental findings. As a mean-field the- 
ory our model can not explain the zero bias anomaly in 
the conductance. However, the parameters of the reso- 
nance state causing the spin-polarization can be used in 
the Kondo model [5]. 
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